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Abstract 

The transmission rate, delivery reliability and network lifetime are three fundamental but conflicting 
design objectives in energy-constrained wireless sensor networks. In this paper, we address the optimal 
rate-reliability-lifetime tradeoff with link capacity constraint, reliability constraint and energy constraint. 
By introducing the weight parameters, we combine the objectives at rate, reliability, and lifetime into a 
single objective to characterize the tradeoff among them. However, the optimization formulation of the 
rate-reliability-reliability tradeoff is neither separable nor convex. Through a series of transformations, a 
separable and convex problem is derived, and an efficient distributed Subgradient Dual Decomposition 
algorithm (SDD) is proposed. Numerical examples confirm its convergence. Also, numerical examples 
investigate the impact of weight parameters on the rate utility, reliability utility and network lifetime, 
which provide a guidance to properly set the value of weight parameters for a desired performance of 
WSNs according to the realistic application's requirements. 
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I. Introduction 

Wireless Sensor Networks (WSNs) consist of a large number of spatially distributed au- 
tonomous sensor nodes with limited computation and sensing capabilities, to monitor physical 
or environmental conditions, and to cooperatively pass their data to a sink. They have been 
extensively applied in many fields, such as battlefield surveillance, environmental monitoring, 
home automation, critical infrastructure protection fll] and so on. 

Recently, there are increasing numbers of network applications, where their performance is 
highly dependent on the high data rate and thus high link capacity requirement. However, the 
link capacity is limited in the WSNs. Thus, many researchers focus on flow/congestion/rate 
control designs to achieve efficient and fair resource allocation in WSNs. The basic framework 
of Network Utility Maximization (NUM) proposed in [2] has been extended to solve flow control 
problem in WSNs flU. Furthermore, the generalized NUM framework proposed in [5] also 
has been used as a tool of cross-layer design in WSNs O. However, all these works mentioned 
above assume that each link provides a fixed-size transmission "pipe" and each user's utility is 
a function of transmission rate only. Furthermore, they don't consider the reliable data delivery 
requirement, and implicitly assume an error-free physical layer, which is hard to achieve in real 
WSNs. 

Since the application performance correlates to the rates of data obtained reliably in WSNs, 
it is vital to guarantee the data delivery reliability requirement in WSNs. There are increasing 
research efforts to improve the reliability: reducing the probability of data loss or error and 
retransmitting data once loss or error occurs. In these work, hop-by-hop recovery [7], end-to- 
end recoveryy (8), and multi-path forwarding [0 are the major approaches to achieve the desired 
reliability. In PSFQ [7], the basic premise is to propagate the segments from source nodes in a 
relatively slow pace and to allow nodes experienced data loss to recover any missing segments 
from immediate neighbors aggressively. In ESRT BE), exploiting the fact that the redundancy 
in sensed data collected in dense WSNs can mitigate channel error and node failure, the sink 
adaptively achieves the expected event reliability by controlling the reporting frequency of the 
source nodes. In RelnForM @, it is proposed to deliver packets at desired reliability by sending 
multiple copies of each packet along multiple paths from sources to sink. Obviously, the data 
transmission rate and the data delivery reliability are two fundamental, yet conflicting, design 
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objectives in WSNs. There is an intrinsic tradeoff between them. However, all these works 
mentioned above didn't consider the intrinsic rate-reliability tradeoff problem. Recent work [fTOl 
has firstly addressed the rate-reliability tradeoff problem explicitly. Through the extended NUM 
framework, where the user utility depends on both transmission rate and delivery reliability, 
the optimal rate-reliability tradeoff can be controlled by adapting channel code rate in each 
link's physical-layer error correction codes. However, it did not take the energy constraint into 
consideration, which is one of the most important constraint in WSNs. 

Typically, sensor nodes are battery-powered, and battery replacement is impossible in many 
sensing applications. Energy is a scarce resource, and WSNs have a finite operational lifetime. 
Hence, network lifetime maximization has been a popular research direction in WSNs, for 
example, ^ has studied the network lifetime maximization problem that jointly considers the 
physical layer, MAC layer and routing layer. [fTT| design an near optimal joint routing-and-sleep- 
scheduling strategy to maximize the network lifetime. [fl3ll propose energy and cross-layer aware 
routing schemes for multichannel access WSNs that account for radio, MAC contention, and 
network constraints, aiming to maximize the network lifetime. [fT4l| propose HYbrid Multi-hop 
routiNg (HYMN) algorithm, which is a hybrid of the flat multi-hop routing and hierarchical multi- 
hop routing, to adequately prolong the lifetime of severely resource-constrained sensor nodes. 
[fT2ll address joint routing and link rate allocation under bandwidth and energy constraints to 
prolong network lifetime and to improve throughput. However, higher data rate leads to greater 
sensing and communication costs across WSNs, resulting in more energy consumption and 
shorter network lifetime. Thus, there is an inherent tradeoff between transmission data rate and 
network lifetime in WSNs. This problem has been extensively studied in recent years [fT5Tl-[fT8Tl. 
but all these works do not consider the reliability requirement in transmitting the data. On the 
other hand, to improve the desired reliability, the approaches, including hop-by-hop recovery 
0, end-to-end recoveryy []8]|, and multi-path forwarding flU, generate more data packets to 
be transmitted, leading to more energy expenditure, and shorter network lifetime. Thus, the 
network lifetime and the data delivery reliability are also two fundamental, yet conflicting, 
design objectives in WSNs. 

It is clear that there is an inherent tradeoff among the data rate, reliability and network 
lifetime: A high data rate can be obtained on a link at the expense of lower delivery reliability, 
which results in more energy consumption and a reduction of network lifetime. Obviously, 
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different applications have very different requirements for rate, reliability, and network lifetime. 
For example, in emergency rescue and disaster relief, it requires a high data rate and reliability, 
but does not have high requirements for the network lifetime. However, in precision agriculture, 
the requirement is to prolong the network lifetime and to improve reliability as much as possible, 
but the data rate is less demanding. Although many works have extensively studied the data rate, 
reliability and network lifetime in recent years separately, so far, no works consider three goals 
together, and study the tradeoff among them. Thus, it is vital to investigate the tradeoff problem 
of data rate, reliability and network lifetime and to design a efficient distributed algorithm to 
achieve the optimal rate-reliability-lifetime tradeoff. 

In this paper, we address the rate-reliability-lifetime tradeoff problem in energy-constrained 
WSNs with link capacity constraint, reliability constraint and energy constraint. First, we in- 
troduce the weight parameters, which combine the objectives at rate, reliability, and lifetime 
into a single objective to characterize the tradeoff among them. However, our new optimization 
formulation for the rate-reliability-reliability tradeoff is neither separable nor convex. It is difficult 
to derive a distributed algorithm that converges to the globally optimal solution. Fortunately, 
through a series of transformations, we convert the formulation into a separable and convex 
optimization problem. Then, the Subgradient Dual Decomposition algorithm (SDD) is applied 
to achieve the optimal solution. Finally, we investigated the impact of different weight parameter 
on the rate utility, reliability utility and network lifetime through numerical examples. 

The rest of this paper is organized as follows. Section II describes the system model. Section 
III formulate the rate-reliability-lifetime tradeoff problem, and transforms it into a separable 
and convex optimization problem. Section IV develops the SDD algorithm to solve the tradeoff 
problem, and proves the convergence of the algorithm. Section V provides numerical examples 
for the proposed algorithms, and illustrates the optimal rate-reliability-lifetime trade-off. Finally 
section VI concludes this paper. 

II. System Model 

In this paper, we consider WSNs consisting of a set of sensor nodes denoted by S = {1, 2..., S} 
and a set of sink nodes denoted by J\f ={1,2,..., TV}. Sensor nodes are battery driven, non- 
rechargeable and irreplaceable. We assume that the sink nodes have enough energy. The sensor 
nodes are the sources that collect data and deliver it to any of the sink nodes, possibly over 



March 18, 2013 



DRAFT 



5 



multiple hops. The WSNs is modeled as a connectivity graph G(V,C) 9 where V = S U A/*, 
includes both the sensor nodes and the sink nodes, £= {1, 2,...,L} represents the set of logical 
links between nodes in the network. We assume that the single-path route is adopted in this 
paper. The key notations used throughout this paper are summarized in Table 1. 



TABLE I 

Summary of key notations 



Notation 


Meaning 


S 


The set of sensor nodes 


C 


The set of logical links 


S{1) 


The set of sensor nodes using link I 


C(s) 


The set of links used by sensor node s 




The incoming links set of sensor node s 


£out(s) 


The outgoing links set of sensor node s 




The set of source nodes that use sensor node 




s as a relay 


St(*) 


The set of sensor nodes that sensor node s 




uses as relays 


](*,*') 


The incoming link of sensor node s on the 




path of sensor node s f 


L out 


The outgoing link of sensor node s on the 




path of sensor node s' 


Is 


The outgoing link that sensor node s uses for 




transmitting its own data 


Q 


The maximum capacity of link I 


™min 
x s 


The minimum data rates for sensor node s 


T max 


The maximum data rates for sensor node s 




The minimum reliability requirement of sensor node s 


j^max 


The maximum reliability requirement of sensor node s 



A. Introducing reliability into NUM framework 

Basic NUM framework assumes that each link provides a fixed size transmission 'pipe' and 
each user's utility is only a function of transmission rate. But in many practical systems adapting 
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the physical layer channel coding or transmission diversity, these assumptions will break down. 
Here, we extended the basic NUM framework to the case that utility of each user depends on 
both transmission rate and delivery reliability, with an intrinsic tradeoff between them. 

Sensor node s sends packets into the encoder of link I at information data rate x s . Packets are 
encoded at the code rate r^ s , where the code rate r^ s is defined by the ratio of the total number 
of useful information bits to the total number of bits transmitted from the encoder per unit time. 
The information bits transmitted from the encoder are sent by the wireless link / at the rate (j)^ SJ 
then (f> liS can be stated as ij5 = x s /r^ s . 

Since the sum of transmission rates of sensor nodes that traverse the link I can not exceed 
the maximum link capacity. Thus, we have 

ses(i) ses(i) 1,8 

The error probability of data transmitted by node s using link / is defined as E(r^ s ), which is 
assumed to be an increasing and convex function of r^ s . 

Let ^ s denote the end-to-end error probability of each node s , then £ s is given by 

£ s = l - J] (l - E(nJ) 

lec(s) 

In general, the error probability of each link is very small, so the end-to-end error probability of 
node s can be approximated as £ s « ^2iec(s) ^( r M)- Let R s denote the reliability of information 
transmitted by sensor node s, then 

* = !-«. ""I-Efcc,.,^) (2) 

Now, we introduce the reliability into the NUM framework. We assume that each sensor 
node s has a utility function U s (x s ,R s ), which is strictly concave increasing functions of the 
information data rate x s and delivery reliability R s . 

B. Network lifetime maximization problem 

In a typical sensor network, sensor nodes have much tighter energy constraints than the sink 
nodes. Hence we will focus only on the energy dissipated in the sensor nodes. Since in most 
types of sensor nodes, communication modules dominate the energy consumption, we ignore 
energy consumed by other tasks such as sensing and data processing. So, we adopt the same 
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simple energy consumption model as in HT8II for the communication module of all nodes. The 
total power dissipation at node s is given by: 

ps = Y Y p7 si Xs ' + Y Y pt *i Xs ' (3) 

iec in (8) s'es(i) lejCout(s) a'eSQ) 

where p r sl is the energy consumption for receiving unit date from link I at sensor node s. In this 
paper, we assume it to be a constant. p l al is the energy consumption for transmitting unit data 
over link I at sensor node s, and is given by: 

Pll = i> + rfal 

where i/j is the electronics energy and c is the amplifier energy, they are constants and depend 
on the function of the physical layer and the environment factors, d s i is the distance of link I at 
sensor node s 9 6 is the path loss factor(2 < 6 < 4). 

We assume that the initial energy of node s is denoted by e s , then the lifetime of node s is 
given by T s = e s /p s . The network lifetime denoted by T is defined as the time at which the 
first node in the network drains out of energy, then T = min 5G< sT s . 

III. Optimal Rate-Reliability-Lifetime Tradeoff 
Now, we formulate the rate-reliability-lifetime tradeoff problem as follows: 

max Y (lsU s (x s , R s ) + (l--f s )wmm seS T s ) (4) 

s 

subject to constraints (pD), ([3]) , and 

R S <1- Y E ( ri ^' seS (5) 

lec(s) 

p s = e s /T s , s eS (6) 
xf n <x s < x^ ax , s e S 
Rf n <R S < i?r x . s e S 
< ri t8 < 1, I e £, 5 e 

7s (0 < 7 S < 1) is a weight parameter to combine different objective functions together into 
a single one. w is a mapping parameter to ensure the objective functions at a same level. In 
this problem, constraint © is reliability constraint for information transmitted by each source. 
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Since the objective function is an increasing function of R s , the inequality constraint © will be 
satisfied with equality at the optimal solution of R s . 

It is very difficult to solve the original problem in a distributed manner since the lifetime 
maximization problem can not be resolved by a distributed manner. We can use a mathematical 
skill to approximate the network lifetime maximization problem for solving it in a distributed 
manner. Considering a general utility function (•) defined by 



V? (x) = 



logx, /? = 1 

■^x 1 -^ p > 1 



1-/3" 

Note that maximizing the minimum rate allocation problem for each source can be approximated 
by maximizing the aggregate utility when the utility function is given in the above form and 
f3 oc [fT9l . The network lifetime maximization problem(i.e., max min sG sT s ) is similar to 
max-min rate allocation problem. We introduce a new utility function Vf (T s ) for each sensor 
node s as a function of its lifetime, which is given by 

vf (T s ) = r L /^ 1 ~ /? 

Then, the maximum network lifetime can be approximated by maximizing the aggregate life 
time utility, i.e., max ^2 aeS Vf (X 5 ). Since the constraint © is not convex and separable, we 
introduce a new variable z s = 1/T 89 which can be interpreted as the normalized power dissipation 
of sensor node s. Then, the constraint © becomes 

p s = e s z s Vs eS (7) 

As a result, the utility function of network lifetime maximization problem has to change corre- 
spondingly: 



Then, the objective function © is transformed into 



Obviously, this objective function is strictly concave. Notice that constraint (OD) in the original 
problem is not satisfied with the properties of separability and convexity, this leads to the original 
problem is neither a convex problem nor a separable one. In the next, we will convert the original 
problem into a separable and convex optimization problem through a series of transformations. 
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First, we introduce a auxiliary variables q jJ5 , which can be interpreted as the allocated trans- 
mission capacity to sensor node s on the link /. Then, the constraint (OD) is decomposed into two 
constraints i.e., 

— <Q >fl ,Ze £,S6S(/) (10) 

V ci s <C h leC (11) 

Notice that the inequality constraint (ITOl) is still inseparable. We take logarithm on both sides 
of this constraint, i.e., logx s — logr/ 5S < logQ j5 . Let x' s = logx s (i.e., x s = e x ' s ), then the 
constraint (ITOl) is changed into 

x' s - logn ?s < logQ jS , l e £,s eS(l) (12) 

Correspondingly, the objective function is transformed into 

W'(x' s , R s , z.) = ls U' s (x' s , R s ) - (1 - % )^ T ^/- 1 

where U r s (x f s) R s ) = U s (e x ' s .Rg). However, U f s (x f s) R s ) may not be a concave function , even 
though U s (x S) R s ) is a concave function. The lemma 2 shown in ifTOll provides a sufficient 
condition for its concavity, under which, the objective function W\x r s) R S) z s ) is also a concave 
function. 

Proposition 1: The rate-reliability-lifetime tradeoff problem is equivalent to the convex prob- 
lem 

max Y, w \ x 'si R s,Zs) 

s 

subject to constraints ( TTT1) . ( IT21) and 

E E P&c^ + E E Plie<' < e s z s ,seS 

iec in (s) s'es(i) ie£out(s) s'es(i) 

I min / / / / max _ o 

x s <x s <x s , s eS 
Rf u <R S < #r x , seS 
< r J>a < l , / e C, se S(l) 
< c hs <Q , I e jC, sg <S(Z) 

where ^ min = logo;™ 11 , and x' s max = logx™ ax . 

Proof First, the constraints © and © can be combined into a single one 



(13) 



= E E ^ + E E p-i 

zec in (s) s'es(z) ze£ „t(s) s'e<s(z) 



March 18, 2013 DRAFT 



10 

By variable substitution x s = e x ' s , the above equality constraint reduces to 

E E P^'" + E E V l sl e<'=e s z s (14) 

Hence, the rate-reliability-lifetime tradeoff problem can be expressed as 

max ^W(x' 8 ,R 8 ,z 8 ) 

s 

subject to constraints ([5]), f fTTl ). ( 1X21 ) and 

E E pS,^ + E E pS,^ = c*. 

iec in (a) s'es(i) iec out (s) s 'es(i) 

< min < x' s < < max , 5 G 5 (15) 

< ri 9a < 1 , / G £, 5 G 
< c M < Q , / G £, 5 G 5(Z) 
Due to the objective of network lifetime maximization (equivalently, minimizing z s 's), it is easily 
known that, at the optimality of the problem (fT3l) . the inequality constraint 

Z Z^ e</ + Z Z v\i^<^zs 

ieCi n (s) s 'es(i) ieC ou t(s) s'es(i) 

must hold with equality. Hence, the problem (fT5l) can be equivalent to (fT3l) which is a convex 
problem. This completes the proof. 

IV. SUBGRADIENT DUAL DECOMPOSITION 

In this section, we will use subgradient dual decomposition approach to solve the problem 
(fT3l) . We write down the Lagrangian function associated with the problem (IT3T) as in the first 
equality of (fT6l) and rearrange the Lagrangian function as in the last equality of Eq. (fT6l) . where 
the new sets S in (s), S t (s) and variables l^ 8 \ l^t \ h are introduced for ease of separation, 
Aj jS , ji s and v s are the Lagrange multipliers which can be respectively interpreted as congestion 
price on link /, reliability price and energy consumption price on sensor node s. 

In Eq. (fT6l) . A s = J2i e c(s) \s > i- e -> the end-to-end congestion price at the sensor node s , and 
p(s',s) _ pr _^_pt ^ j e ^ ^ p 0wer dissipation of relaying unit data from node s at node 
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£ (x', R, r, c, z, A, fi, v) 

= ^W'(x' s ,R s ,z s ) + Y: E A,, s (logq, s + logr,, s -<) + E^ (l- E E(r ltS )-R s 
s i seS(i) s \ lec(s) 

+ Ev s [e s z s - E E P r si^' ~ E E P**e< 
* \ iec in (s) s'es(i) ieC ou t(s) s 'eS(i) 



^\W(x' s ,R s ,z s )-x' s E Ks ~ VsRs + v s e s z s - E "sP r j s ,s>)e x ' s ' ~ E 



(s) 



+ E<i E (A/, s (logc M + logr i>s ) - fj, 8 E(ri )S )) } + 
se5(«) 



: E S r s, z s) ~ AV S - n s R s + i/ a e a z a - E v s'f , .(.',.> - E v s'V t , , (s ', s) - ^hli, I 



+ ES E (A*, s (logq ;S + logr i;S ) - fjL s E(rt, s )) > +E/ i * 
= E \ W'(x' s , R s , z s ) - X s x' s - fi s R s + v s e s z s - e< E *V ( P r , ,(.',.) + P*, ,(.',.) ) 

+ ES E (A/, s (logc /jS + logr iiS ) - fi s E(r l>s )) \ + E/^ 
= E S w '( x 'si R s, z a ) ~ AV S - fi s R s + ^ s e s 2; s - e x '° E ^'P (s '' s) - v s e x '*ip\ ls 

s [ s'eSt(s) 

+ ES E (Aj, s (logq ;S + logr i;S ) - fi s E(r l>s )) \ + Ep* 
z I se5(0 s 



v s e< V l sl 



(16) 



s' ,s ; G <St(s) , s G <S . The Lagrange dual function is given by 

G (A, i/) = max £ (x , R, r, c, z, A, /x, i/) 
subject to x' min < x' r< x' max 

R min ^ R ^ R max (17) 

■< r ^ 1 

c € C 

where C = {(c,, a )Z € £, s € 5(Z)| E s6 s ( z) C M < Q, Z G £, < c M < C,, Z G C, s G <S(Z)}. 
The dual problem corresponding to problem ([13) is then given by 

min G(X,u,u) (18) 
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The key of the dual-decomposition algorithm is solving the problem (ITTI) separately and 
distributively. With the separation in Eq. (fT6l) . maximization the Lagrangian over (x, R, r, c, z) 
can be done in parallel at each sensor node s 

max W'(x f 8 , R s , z s ) — \ s x' s — fi s R s 

+v s e s z s - e x ' s Y ^s'P^'' s) ~ Vs^'s^si 

s'eMs) (19) 

g -j- jJ min <^ x max 



RT n <R H <R 



max 



and at each link Z 



max E (A/, s (logQ ?s + logr M ) - n s E(r hs )) 
ses{i) 

s.t. E % a <d 

seS(l) UU) 
< Q, s < Ci , s G <S(Z) 

< r M < 1, s e 

The problem (l20l) can be further decomposed into two sub-problems as follows: 
Link-layer sub-problem 

max E \ 5 logQ j5 

s65(Z) 

s.t- E Q >fl <Ci (21) 

8ES(l) 

< q, s < Q , 5 G <S(Z) 
and physical-layer sub-problem for sensor node 5, 5 G <S(Z) 

max A M log r M - fi s E(r^ s ) 
s.t. < r M < 1,8 G <S(Z) 
Once the problem ([T71) is solved, the subgradients of the dual function with respect to the dual 
variables can be calculated easily and the dual variables for solving the dual problem (IT81) can 
be iteratively updated by using subgradient projection method [20J as follows: 
Congestion price update at each link Z, Z G C, s G <S(Z) 

A M (t + 1) 

= [\ s (t) - 5 (t) (logc M (t) + logr M (t) - + (23) 

= [A M (t) - 5 (t) (log Q , s (t) + logr M (t) - logx s (t))] + 
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Vs (t + 1) 



v s (t) -0(t)[ e s z s (t) -EE f sl e<> {t) "EE Plie<' {t) 

ieCi n (s) s 'eS(i) ieC ou t(s) s 'eS(i) 



v s (t)-$(t) \e s z s (t)- J2 e 



</(*) 



s'eSin(s) 



p 



As,s>) +P t .(s,s') 
,L in bi out 



pX s (t) j^t 



is s (t) - $(t) ( e s z s (t) - £ x s ,(t)p^ - x 8 (t)p t 

s'es in (s) 



(25) 



Reliability price update at each sensor node s, s e S 

(t + 1) 



fi s (t)-((t) l- E E(n 9a (t))-R a (t) 

\ lec(s) 
[fi s (t)-((t) (R s (t)-R s (t))] + 



(24) 



where R s (t) 
node s. 



E E( r i,s(t)) with an interpretation of end-to-end reliability at a sensor 



Z<E£(s) 



The energy consumption price v s is updated according to (1251) (see the top of the next page), 
where p( s > s "> = p r (s s/) + p l (s s/) , i.e., the power dissipation of node s for relaying unit 
data from node s' , s' e S in (s) , 5 G 5 . In the above formulas, [w] + = max{0, w} , and 5 (t) , 
( (t) and i? (t) are positive scalar step size. 

We will summarize the distributed algorithm for rate-reliability-lifetime tradeoff as follows, 
where each sensor node and each link will solve their own problems with only local information. 
The information exchange between sensor node and link in the distributed algorithm SDD is 
given in Fig. 1. 



SDD: Subgradient Dual Decomposition Algorithm 



at each iteration t 

at each sensor node s 

1) Rate, reliability and lifetime update: 
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• Node s receives the link congestion price \i l8 (t) of link I from the network, I e C(s) 9 
and calculates A s (t) according to A s (t) = $^ze£(s) \s(*)- 

• Node s receives the energy consumption price v s > (t) of sensor nodes s' who relay data 
packets for node s, s f e S t (s). 

• Node s locally solves the problem (IT91) for the given /i s (£), and A s (t). Up- 
date the information rate x s (t + 1) (where x s (t + l) = e x ' a ^^\ information reliability 
R a (t + 1) and node lifetime T s (t + l)(where T s (t + 1) = l/z s (t + 1) ). 

• Broadcasts the new rate x s (t+ I) to links that sensor node s uses. 

2) Reliability price update: 

• Node s receives code rate r^ s (t) of link I, I 6 £(s), and computes the end-to-end 
reliability for given R s (t) = 1 — Y^iec^s) E{r^ s (t)), then, updates its reliability price 
according to Eq. d24l) . 

• Broadcasts the new reliability price ji s [t + 1) to the links that sensor node s uses. 

3) Energy consumption price update: 

• Node s receives x 8 t(t) of the sensor nodes that use node s relaying their data packets, 
s' G Si n (s), and updates its energy consumption price according to Eq. d25l) . 

• Broadcasts the new energy consumption price v s (t + 1) to the nodes that use node s 
relaying their data packets. 

at each link /: 

1) Auxiliary variables q jS update: 

Link / update Q ?s (t) by locally solving the link-layer problem (I2T1) for given A/ j5 (i). 

2) Code rate r^ 5 update: 

• Link I receives reliability price n 8 {t) of sensor node s, s e <S(Z) , then, update 

+ by locally solving the physical-layer problem d22l) for given /i s (t) and 
A, )fl (t). 

• Broadcasts new code rate r^ s (t + 1) to the sensor nodes that use link /, s e 

3) Congestion price update: 

• Link I receives data at rate x s (t) from the sensor nodes that use link /, s e S(l) , and 
updates its congestion price according to Eq. (1231) . 

• Broadcasts the new congestion price A/ j5 (t + 1) to the nodes that use link /, s e <5(Z). 
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(I2TT) shows the global convergence of the proposed dual-decomposition based algorithm. For 
convenience, we state the convergence result in Theorem 1. 

Theorem 1: By the above distributed algorithm, dual variables (\(t), fj,(t) 9 v(t)) converge 
to the optimal dual solutions (A*, i/* ), if the setpsizes are chosen such that 5 (t) — » 

CO CO CO 

0, E s (*) = 00 > C (t) -> E C (*) = oo and i? (t) ^ 0, £ # (t) = oc. 
t=i t=i t=i 

On the other hand, since the tradeoff problem (IT3T) is a convex optimization problem, by 
Lagrange-duality theory [|201 . we conclude that, the corresponding primal variables (x*, R*, z*, c*, r* 
) are the optimal solutions of the problem (fT3l) . 

V. Numerical Examples 

In this section, we present numerical examples for the proposed algorithm by considering a 
sensor network, shown in Fig. El In this network, we have six sensor nodes indexed 1-6, and one 
sink node. They communicate through seven links. The locations of sensor nodes are randomly 
generated over a 100m x 100m square area. 

We set the objective function W(x SJ R S) z s ) in the following form 

^.l — ct _^min(l — a) 

W[X S) R S) Z s ) 7s^5 max(l-a)_^ min(l-a) 

X s X s 

, ^ n _ m \ Ri- a -Rf n(1 - a) 

~T IS\*- Sks/ D max(l-«) D min(l-a) 
H s —rC s 

-(i-Ts)^/" 1 «>0, a^l 
where <p a (0 < <p a < 1 ) is the weight parameter that weights the relative importance of 
information data rate and reliability. The three parts of the expression of the objective function 
represent the rate utility, reliability utility and network lifetime, respectively. E{r^ s ) is assumed 
to be of the following form: 

E ( r i,s) =\ ex P (-« (i - n, s )) 

where k is the code block length used by the encoder. In our experiments, the constant parameters 
are set as follows: xf n = 0.1 Mbps, x™ ax = 2.0 Mbps, Rf n = 0.9, i?™ ax = 1, a=l.l, (3=9 
and the mapping parameter tu=3.2768 x 10 32 . 

We assume that there is a routing mechanism in place to find a route for each sensor node. 
The routes of each flow have been drawn in Fig. |2l The capacity of links a — g are set to 
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be C\ = [2; 3; 2.5; 3; 4; 2.5; 4] (Mbps). The initial energy of sensor nodes are set to be e = 
[3000; 2800; 2500; 2200; 2600; 2000] (J). The power dissipation at node s is determined by @ 
in section II, where ij) = 50nJ/b , s = 0.0013pJ/6/ra 4 , 6=4 and p r sl =50nJ/b [|22t 

A. Convergence performance of SDD 

In the section, we depict the convergence performance of SDD algorithm. We only show the 
convergence results of nodes rates, the convergence results of others are omitted due to space 
limitations. Both the weight parameter ^ s and tp s are set to 0.8, and the rates of sensor nodes 
are shown in Fig. [3j Notice that the rates of all sensor nodes change sharply at the beginning 
of the iteration, and then converge to the optimal solution after about 200 iterations. 

FigJH depicts the convergent behavior of total utility solved by SDD. The blue dotted line 
denotes the optimal value of the total utility solved by the centralized algorithm. Obviously, the 
values of the total utility solved by SDD converge to the optimal value after about 200 iterations. 

Then, both the weight parameter ^ s and (p s are increased to 0.97, in which the reliability and 
network lifetime is almost out of consideration. The convergent results of nodes rates and total 
utility of SDD are shown in Fig. [5] and Fig. [6l respectively. Notice that the optimal rates of all 
the sensor nodes are larger than those in Fig. |3] since we increase 7 S and cp s to obtain more rate 
utility while reducing the data reliability and network lifetime. 

Furthermore, both the weight parameter j s and (p s are set to 0.5 , in which the network lifetime 
is more consideration than the rate and reliability. The convergent results of nodes rates of SDD 
are shown in Fig. [71 Compare with Fig. |3] and Fig. |5] the optimal rates of all the sensor nodes are 
very small. In this case, the sensor nodes can save more energy to extend the network lifetime. 
FigE] depicts the convergent behavior of total utility. Compare with Fig. |4] and Fig. |6l We can 
see that the total utility has a significant reduced. 

B. The impact of weight parameters on the rate utility, reliability utility and network lifetime 

In this section, we will investigate the impact of weight parameters j s and (p s on the rate 
utility, reliability utility and network lifetime. The inherent tradeoff between rate utility and 
reliability utility can be observed from Fig. where j s =l (i.e., network lifetime isn't taken 
into account), and (p s ranges from to 1. We can see that as (p s increases from to 1, the rate 
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utility increases, but the reliability utility decreases; and the larger the (p s , the larger drop in 
reliability utility. 

Fig. [10] shows the inherent tradeoff between rate utility and network lifetime, where (p s =l 
(i.e., reliability isn't taken into account), and j s ranges from to 1. The network lifetime is 
shown in seconds. From Fig. [lOl we observe that as 7 S increases from to 1, the rate utility 
increases, while the network lifetime decreases and drops sharply when j s is small. So, there is 
a evidently tradeoff between rate utility and network lifetime. 

In the realistic applications, according to the actual requirements, we properly determine the 
value of weight parameters j s and (p s to achieve a desired performance at the transmission rate, 
delivery reliability and network lifetime. 

VI. Conclusions 

We have studied the rate-reliability-lifetime tradeoff in WSNs with link capacity constraint, 
reliability constraint and energy constraint. Our study is the first to jointly consider rate, relia- 
bility, and network lifetime in this rigorous tradeoff framework. Our optimization formulation 
for the rate-reliability-reliability tradeoff is neither separable nor convex. We convert the new 
formulation into a separable and convex optimization problem through a series of transformations. 
A distributed SDD is developed for solving the rate-reliability-lifetime tradeoff problem. We 
also demonstrate that the convergence speed of SDD through numerical examples. Finally, We 
investigated the impact of different weight parameters j s and (p s on the rate utility, reliability 
utility and network lifetime through numerical examples. We can select the appropriate value 
of weight parameters according to the actual requirements to achieve a desired performance of 
WSNs. 
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Fig. 1. Diagram for the distributed algorithm. 




Fig. 2. Topology of wireless sensor network. 
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Fig. 3. The convergent performance of node rates in SDD with 7 S =0.8 and (p s =0.S . 
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Fig. 5. The convergent performance of node rates in SDD with 7 S =0.97 and (p s —0.97 . 
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Fig. 7. The convergent performance of node rates in SDD with 7 S =0.5 and (p s —0.5 . 
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